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The magnetic field and geometry of the oblique shock in 
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ABSTRACT 

We investigate the brightest regions of the kpc-scale jet in the powerful radio 
galaxy 3C 346, using new optical Hubble Space Telescope (HST) ACS/F606W 
polarimetry together with Chandra X-ray data and 14.9 GHz and 22.5 GHz 
VLA radio polarimetry. The jet shows a close correspondence in optical and 
radio morphology, while the X-ray emission shows an 0.80 ± 0.17 kpc offset from 
the optical and radio peak positions. Optical and radio polarimetry show the 
same apparent magnetic field position angle and fractional polarization at the 
brightest knot, where the jet undergoes a large kink of almost 70 degrees in 
the optical and radio images. The apparent field direction here is well-aligned 
with the new jet direction, as predicted by earlier work that suggested the kink 
was the result of an oblique shock. We have explored models of the polarization 
from oblique shocks to understand the geometry of the 3C 346 jet, and find that 
the upstream flow is likely to be highly relativistic {(3 U — 0.91~^ q 7 ), where the 
plane of the shock front is inclined at an angle offy = 51±ll degrees to the 
upstream flow which is at an angle 8 = 14+y degrees to our line of sight. The 
actual deflection angle of the jet in this case is only 22 degrees. 

Key words: galaxies: jets - galaxies: active - galaxies: individual(3C 346) - 
magnetic fields - polarization 



1 INTRODUCTION 

Radio observati ons of kpc-scale jets in active galaxies 
have shown fe.g. lBridle fe Pcrlcv 1983) that the emissions 
tend to be highly polarized and exhibit a characteristic 
synchrotron spectrum. Observations of nearby sources at 
shorter wavelengths often show a close correspondence 
to the radio morphology and characteristic spectral- 
energy distributions, suggesting that synchrotron radia- 
tion from energetic particles also dominates the emission 
in th e optical (e.g. [Keel 1988: iFedorenko fc Courvoisiei] 



1996) and X-rav |Worrall. Birkinshaw 



2001 



Hardcastld 

[Hardcastle. Birkinshaw fc Worralll 1200 ll ). High- 
energy electrons (7 = E/mc 2 ~ 10 7 - 10 8 ) radiate 
at optical and X-ray wavelengths and have very short 
synchrotron-emitting lifetimes (10 2 - 10 3 years in a mag- 
netic field of a few nT), and must therefore be reacceler- 
ated locally to maintain the observed emissions. Lower- 
energy electrons have much longer lifetimes (~ 2 x 10 6 
years, for 7 = 10 4 and B = 3 nT) and are detected pri- 
marily at radio wavelengths. 

Polarization studies of the synchrotron emission from 



kpc-scale jets provide valuable information about their 
structures and magnetic configurations: the observed po- 
larization direction gives an emission-weighted measure 
of the apparent magnetic field direction along the line of 
sight, while the fractional polarization indicates the or- 
dering of the field. Multi-band imaging and polarimetry 
are therefore powerful tools in the study of jets, allowing 
us to use the emission across the electromagnetic spec- 
trum to probe populations of particles at different ener- 
gies. 

The powerful FR II radio galaxy 3C 346 lies at 
redshift z — 0.16 and has a 1 7 th -magnitude elliptical 
host at optical wavelengths (e.g. lLaing. Riley fc Longaiij 
198 3]). Althoug h origi nally classified as an FR I galaxy 



by lLaing et al.l (|l983l l. radio maps of the source show 
a highly luminous core surrounded by diffuse lobe emis- 
sion and a highl y-distorted eastern jet exhibiting a brigh t 
knot structure JSpencer et al.lll99ll : ICotton et alJll995t ): 
ISpencer et al. I i|l99ll ) argued that these features could 
be best explained by an FR II radio galaxy being fore- 
shortened by rel ativistic beam i ng at a small angle to the 
line-of-sight, and lCotton et all (|l995l l used the radio-core 
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dominance and jet/counter-jet ratio to deduce an angle 
to the line-of-sight 8 < 32 degrees and jet speed /3 > 0.8. 
The jet was first detected i n the optical using g round - 
based observations made bv iDev fe van Breueell (1 19941 ). 
Subsequent Hubble Space Telescope (HST) observations 
showed striking simil arities between the optical and ra- 
dio j et morphologies ijde Kofi" et al ] ll996t Ide Vries et al] 

Il997h . 

IWorrall fc Birkinshawl (|2005h reported the first de- 
tection of X-ray radiation from the brightest knot in the 
jet of 3C 346 using high-resolution images from the Chan- 
dra X-ray Observatory, finding that the peak of the X- 
ray emission does not exactly match the peak of the ra- 
dio and optic al emission. The availa b le mu ltiwavelength 
data allowed IWorrall fc Birkinshawl l|2005l ) to conclude 
that the X-ray emission process is synchrotron and to dis- 
cuss models for the jet. They suggested that the bright 
knot and accompanying sharp kink in the jet were the 
result of an oblique shock, and predicted that this should 
generate a relatively simple polarization pattern with the 
magnetic field direction at position angle ~ 20 degrees. 
This is consistent with the high (17 per cent) radi o po- 
larization measured bv lAkuior fc Garringtonl l|l995l ). but 
sufficiently good high resolution optical and radio po- 
larimetry were not available to probe the magnetic field 
structure until now. The polarization data allow us to test 
the prediction that the magnetic field vectors should be 
largely aligned with the direction of the putative oblique 
shock. 

In this paper we re-exami ne the 3C 346 j et, u sing 
new HST optical polarimetry (IPerlman et al] I2006T ) in 
conjunction with previously unpublished 14.9 GHz and 
22.5 GHz VLA radio polarimetry, and the Chandra X- 
ray data. This work probes the jet in sufficient detail 
for us to test the predictions of Worrall fc Birkinshawl 
(2005) and distinguish between different shock models in 
3C 3 46. We adopt mod els for oblique and conical shocks 
from lCawthornei l|2006h . and use these to explore the ge- 
ometry of the shock and the upstream flow speed. Sec- 
tion [2] describes the multi-wavelength observations we 
used. We present our results in Section [3] and in Sec- 
tion |3] we examine models of the magnetic field at the 
shock. Throughout, we adopt values for the cosmologi- 
cal parameters Ho = 70 km s _1 Mpc -1 , Qao ~ 0.7 and 
fimo = 0.3. Spectral indices are defined as the slope of the 
power-law relating the flux S„ at a frequency v according 
to S v cx v~ a . At the redshift of 3C 346, one arcsecond 
corresponds to a projected distance of 2.76 kpc. 



2 OBSERVATIONS 
2.1 Optical data 

All observation details are summarised in Table [T] Opti- 
cal HST images of 3C 346 were obtained on 2003 August 
19 with the ACS/WFC instrument using the wide-band 
F606W filter and the three optical polarizers, each in 
a 502 second exposur e. Optical polar i metry on 3C 346 
was first presented bv lPerlman et all (|2006l ). Data were 
flat-fielded and bias-corrected using standard packages in 
iraf/stsdas, and CRREJECT was run to remove cosmic- 



ray events (see lPerlman et al"1l2006l . 1 19991 . for further de- 
tails ). The ACS imag es were combined using multidriz- 
zle (|Koekemoer et al.ll2002l), allowing for corrections in 
the chip geometry (| Anderson fc Kinell2004l ). Stokes /, Q 
and U images were produced and combi ned following the 
guide lines in the ACS Data Handbook l|Pavlovskv et al] 
120051 ). The optical percentage-polarization and apparent 
Magnetic Field Position Angle (MFPA) images were com- 
puted using 100 x (Q 2 + U 2 )i/I and § tan" 1 (U/Q) + 90 
degrees, respectively. Corrections were needed for the op- 
tical MFPA to include the HST roll-angle (285.9 de- 
grees for this observation) and the offset given by the 
camera geometry (— 38.2 degrees for the WFC, from 
iPavlovskv et al.H2005h. The Ricia n bias in the percent- 
age polarization 1 SerkowskHll962l ) was corrected using a 
script adapted from the Space Telesc ope European Co - 
ordinating Facilit y (ST-ECF) package dHook et al.ll20odh 
in iraf, following IWardle fc Kronberg (|l974l ). The script 
uses a 'most-probable value' estimator, excluding pixels 
where the signal-to-noise ratio is < 0.5, or where the most 
probable value in the percentage polarization is negative 
or > 100 per cent. Small corrections were made for the 
parallel and perpendicular transmittance of each filter. 

Galaxy subtraction was performed using the 
tasks ellipse and bmodel in iraf, and im- 
CALC in STSDAS. The HST point-spread function 
(PSF) was modelled for the ACS/F606W instru- 
ment/filter combination using tinytim version 6.3 
( |http: / / www.stsci.edu/software/tln ytim/tinytim. html ) , 
and measured to have a FWHM of 0.156 arcsec. The 
optical images were subsequently read into AIPS for 
analysis. 



2.2 X-ray data 

The ACIS-S instrument on Chandra was used to ob- 
serve 3C 346 on 2002 August 3 (ObsID 3129) in 
FAINT imaging mode, using chip S3 with a 0.84-second 
frame time and a 256-row subarray. The level- 1 X- 
ray data were reprocessed to take advantage of more 
appropriate calibration products. New bad-pixel files 
were created following the analysis 'threads' (online at 
http://cxc.harvard.edu/ciao/threads/), and a new level- 
1 event file was created with pixel-randomisation re- 
moved. The event file was filtered for good grades with 
values 0, 2, 3, 4 and 6, and time intervals were excluded 
where the background rate was larger than 3-sigma from 
the mean. The mean background rate was 0.3 counts per 
second, and the exclusion threshold was 0.4 counts per 
second. The new level-2 event file contained 41129 s of 
good data, selected from the original 44765 s observa- 
tion. 

The X-ray image was resampled to the same (0.05 
arcsec) pixel scale as the optical data, which is much 
smaller than the native (0.492 arcsec) pixel scale or spa- 
tial resolution of Chandra. The resampled image was then 
smoothed using a two-dimensional circular Gaussian with 
FWHM 0.35 arcsec for the figures in this paper. The re- 
sulting image resolution of ~ 0.7 arcsec FWHM is domi- 
nated by the width of the PSF. 
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Figure 1. Top Panel: 22.5 GHz radio intensity contours and polarization vectors (rotated through 90 degrees), overlaid on 
a heavily clipped greyscale Chandra X-ray image smoothed with a 0.35 arcsec (FWHM) circular Gaussian. Contour levels 
are 0.36 X (1.2,4,8,16,32) mjy beam -1 , and a vector of length 0.15 arcsec corresponds to 100 per cent polarization. 
The angular resolution of the radio image is 0.118 X 0.108 arcsec. Bottom Panel: Our highest-resolution radio map, created 
using a restoring beam of 0.092 X 0.080 arcsec. The figure shows the region around the brightest jet knot superimposed on 
a colour scale of the X-ray emission. The offset between the X-ray and radio peak positions is clear. 



2.3 Radio data 

In order to obtain a high resolution radio map suitable for 
comparison with the HST images, we reprocessed archival 
K-band (22.5 GHz) and U-band (14.9 GHz) VLA data of 
3C 346 using the MPS package from NRAO. A total of 
91.5 minutes of data were obtained at each frequency in 
the A-configuration over two days, starting 1996 Decem- 
ber 14, and the VLA was used in fast-switching mode at 



22.5 GHz to keep track of the rapid phase variations that 
often affect high-frequency interferometer data. 

The data were read into aips and bad scans were 
flagged and removed. Baseline corrections were avail- 
able from the NRAO logs, so these were applied using 
the task CLCOR. Flux-density and polarization calibra- 
tion used 3C 286, and 1658+076 was used for phase 
referencing during both observations. The calibration 
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Table 1. Details of the radio, optical and X-ray observations used in this study of the jet in 3C 346. Integration times are 
given for the total observation period, before processing. 



Instrument, Band/Filter 


ID 


Integration Time (s) 


Date 


"min (Hz) 


t'max (Hz) 


VLA A-config, U-band 


A0127 


5.49 X 10 3 


1996 Dec 14 


1.49 x 10 10 


1.50 x 10 10 


VLA A-config, K-band 


A0127 


5.49 x 10 3 


1996 Dec 14 


2.24 x 10 10 


2.25 x 10 10 


HST ACS/WFC1, F606W - 


- POL 9847 


1.51 x 10 3 


2003 Aug 19 


4.04 x 10 14 


6.76 x 10 14 


Chandra ACIS, 0.3-6 keV 


3129 


44.77 x 10 3 


2002 Aug 3 


7.24 x 10 16 


1.45 x 10 18 



procedure was performed in two passes, following the 
guidelines for reducing high-frequency VLA data (in ap- 
pendix D of the AIPS Cookbook, available online at 



http://www.aips.nrao.edu/cook.html). Phase-only cor- 



rections were applied in the first pass to minimize decor- 
relation in the calibrator scans, and only then was the 
full amplitude-and-phase calibration applied in the sec- 
ond pass to determine the flux density scale. The primary 
calibrator (3C 286) is resolved by the VLA at high fre- 
quencies and in extended configurations, so both passes 
used a clean-component model of 3C 286 to represent 
the small-scale structure on the longer baselines. A good 
calibration was obtained by tracking the rapid phase vari- 
ations using a sufficiently short (3 second) interpolation 
time in the calibration table for the K-band data. The 
bright core of 3C 346 also provides an excellent anchor 
for self-calibration, allowing the data to be phase and 
amplitude self-calibrated to good convergence in Stokes 
I. Stokes Q and U images were then produced from this 
calibrated data set. 

All the radio images were made on a 0.02 arcsec grid. 
Low-noise images of 3C 346 were produced using the K- 
band data with near-natural weighting (robust = 4), to 
give an elliptical Gaussian beam of size 0.118 x 0.108 arc- 
sec and major axis in position angle 90 degrees. This gives 
a noise level in Stokes L of ~ 75 /iJy beam -1 . Our highest- 
resolution K-band image has a beam of 0.092 x 0.080 
arcsec. The U-band data were imaged with a larger 
(0.130 x 0.120 arcsec) beam, and have a lower noise level 
of ~ 40 pJy beam" 1 . The percentage-polarization and 
apparent MFPA images were produced as described in 
section 12.11 The Ricia n bias in the percentage polariza- 
tion (|Serkowskil [19621 ) was corrected by measuring the 
noise levels in the Stokes Q and U images and using the 
POLC algorithm in the AIPS task COMB. Component fluxes 
were calculated by fitting two-dimensional Gaussians us- 
ing the task IMFIT in MPS. 

A radio spectral-index map was produced by com- 
bining the 14.9 GHz and 22.5 GHz VLA data. Both data 
sets were matched in resolution and carefully aligned be- 
fore combining them. The 14.9 GHz data was shifted by 
(-0.009, 0.016) arcsec relative to the 22.5 GHz data to 
ensure good alignment of the core in each image. 

To compare the radio and optical data, the radio im- 
ages were first regridded to match the 0.05 arcsec pixel 
scale of the optical images, and then convolved with an el- 
liptical Gaussian to give the radio beam the same FWHM 
(0.156 arcsec) as the optical PSF. 

Since the absolute astrometry of radio images from 
the VLA is more accurate than either the optical or 
X-ray data, all the cores were registered to the coordi- 



nates of the radio peak (at RA = 16 h 43 m 48:600, Dec. = 
+ 17°15'49'.'44 in J2000) prior to any analysis. The opti- 
cal images were shifted by (—0.17, —0.17) arcsec and the 
X-ray by (0.05,0.06) arcsec in (RA, Dec). 



3 RESULTS 

3.1 Jet morphology 

The eastern jet of 3C 346 shows a variety of knots 
and a dramatic kink in the radio and optical images, 
shown in Figures Q] - [2] The knots are labelled from 
A- G with increasi n g dist ance from the nucleus, follow- 
ing [Perlman^et^al] (2006). Knots A1-A3 are very faint: 
knot Ai is not detected in our 22.5 GHz radio image, but 
we do detect it us i ng th e 14.9 GHz data. 

ICotton et al.l l| 19951 ) used radio Very Long Baseline 
Interferometry (VLBI) data to show that the jet is ini- 
tially straight, from the point at which it is visible against 
the bright galactic-scale emission until it curves to form 
the brighter knots B and C. At knot C the jet undergoes 
an apparent kink of ~ 70 degrees, although this does not 
disrupt it: the jet continues to curve and sh ows additional 
knots as i t moves into the eas tern lobe. Ide Koff et alj 
ljl996h and Ide Vries et all l| 19971 ) noted the striking simi- 
larity between the radio and optical jet morphologies, and 
our matched HST and VLA images of the jet are also in 
excellent agreement (Figure [5} . However, we note that 
the jet appears to be somewhat more 'knotty' in the op- 
tical than in the radio: the optical data show that knot C 
has a more extended southern edge, and all the knots in 
the optical image appear to be more distinct than in the 
smoother radio jet. This phenome non has been observ ed 
in other kiloparsec-scale jets (e.g. iBahcall et al.lll995l . in 
the case of 3C 273). 

IWorrall fc Birkinshawl (|2005l ) reported the discovery 
of X-ray synchrotron emission associated with knot C. 
With the improved X-ray calibration and radio maps that 
we are using, we confirm their claim of a positional off- 
set between the X-ray and radio peaks. The X-ray im- 
age shows that the centroid of X-ray peak, labelled Xbc 
in Figured] is located at RA = 16 h 43 m 48!731, Dec. = 
-r-17 o 15'49 / .'80 (measured using the zhtools software), 
which is 0.29 arcsec away from the peak of the radio emis- 
sion (at RA = 16 h 43 m 48!750, Dec. = +17°15 , 49 , .'78). 
With ~ 70 X-ray counts in this region, the peak should 
be located to better than 60 mas using Chandra, so the 
offset in peak positions is significant at the 4cr level. 

The integrated 22.5 GHz radio flux density at knot C 
was measured to be 54.9 ±0.1 mjy. Its deconvolved size 
is 0.11 x 0.07 arcsec 2 , corresponding to an ellipsoid with 
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Table 2. Polarimetry of the knots in the jet of 3C 346 at 
optical and radio wavelengths, showing the percentage polar- 
ization at 1.3 cm (P ra d)i 606 nm (P op t), and the difference in 
apparent magnetic-field position angle (optical-radio), A x , in 
degrees. 



Knot Prad ' Popt ' A x , degrees 

per cent per cent 



Ai 








A 2 


65 ± 10 


40 ±8 


25 ± 10 


A 3 


65 ± 10 






B 


25 ± 6 


18 ±8 


12 ±7 


C 


20 ±4 


27 ±7 


2±2 


D 


14 ±7 


20 ±6 


-24 ± 11 


E 


30 ± 10 







major and minor axes 300 p c and 200 pc, respectively. 
IWorrall fc Birkinshawl (|2005h measured a 1 keV flux den- 
sity at knot C (Xbc in Figure [TJ of 1.6 ± 0.2 njy, and 
found the X-ray spectrum consistent with an unabsorbed 
simple power law of spectral index a = 1.0 ± 0.3 (for a 
spectrum S v oc v~ a ). The radio spectral index variations 
were mapped using the 14.9 GHz and 22.5 GHz VLA 
data. The core shows a spectral index of a = —0.2 ± 0.1, 
consistent with self-absorbed emission, and the radio 
spectral index around knot C shows low-significance vari- 
ations from a — 0.4 to a = 0.6. 

3.2 Jet polarimetry 

The structure of the jet's apparent magnetic field is in- 
dicated in the optical and radio maps in Figure [2] where 
the polarization vectors have been rotated through 90 de- 
grees. Although this is common practice, these vectors do 
not necessarily show the true projected magnetic field po- 
sition angle (M FPA), which may shift due to relativ istic 
aberration (e.g. iLvutikov. Pariev fc Blandfordlliool ). 

Error maps of the radio and optical MFPA were pro- 
duced by measuring the noise on the Stokes Q and U im- 
ages. In addition to random noise, the optical ACS data 
suffer from a 3-degree system atic uncertainty in the po- 
larization position angle (see iPavlovskv et al. 1 120051 . for 
further details). Polarization vectors are plotted in Fig- 
ure [2] only where the statistical errors in the apparent 
MFPA are less than 10 degrees: regions where the errors 
are larger than 10 degrees are shown without vectors and 
are ignored in this analysis. It should be noted that statis- 
tical errors in the apparent MFPA are often much smaller 
than this, and only about 2 degrees in the brightest knot. 
The differences in polarization position angle between the 
14.9 GHz and 22.5 GHz radio data are small (generally in 
the range — 4 degrees) and, within errors, are consistent 
with zero. 

Figure [5^, (colour-scale) shows the ratio of the per- 
centage polarization detected in the radio and optical 
bands, where values near zero indicate relatively little po- 
larized radio flux. The data show a weak trend towards 
increasing optical/radio polarization ratio with distance 
along the jet: out to knot C, the jet is relatively more po- 
larized in the radio than in the optical. Near the base, in 
knots A2 and A3, the radio polarization is high (65 ± 10 



per cent) but this falls to 25 ±6 percent in knot B. In the 
optical, knot A2 and knot B are 40 ± 8 per cent polarized 
and 18 ±8 per cent polarized, respectively. No significant 
optical polarization is detected in Ai or A3. The radio 
polarization continues to fall, reaching 20 ± 4 per cent 
at knot C and 14 ± 7 per cent at knot D. In contrast, 
the optical polarization does not change much: 27 ± 7 per 
cent at C, and 20 ± 6 per cent at D. Knots E and F are 
weakly detected in the optical but show no significant po- 
larization, and little optical jet emission is detected after 
knot F. 

The radio emission and polarization structure con- 
tinue after knot D, however. We note in particular the 
asymmetric polarization seen across knots D and E, 
which could indicate a helical structure in the magnetic 
field. The north-western side of knot E is highly polarized 
in the radio (30 ± 10 per cent), but the south-eastern side 
of knot E shows no significant polarization. The radio 
polarization is weaker at knot D (as noted above), but it 
too shows a similar asymmetry across the jet. 

The colour-scale in Figure [2jb) shows differences (in 
degrees) between the apparent MFPA seen in the optical 
and radio data. It is clear that the apparent MFPA differ- 
ences are small throughout the jet and, as shown in Table 
2, the optical and radio polarization directions are partic- 
ularly well aligned in the bulk of knot C. The radio and 
optical MFPAs in the knot A complex are broadly paral- 
lel to the jet, although they are misaligned by 25 ± 10 
degrees in knot A2. In knot B the apparent MFPA re- 
mains almost unchanged in the radio, while the apparent 
optical MFPA rotates through a larger angle to reduce 
the misalignment: the differences in knot B are less signif- 
icant, with a misalignment of 12 ± 7 degrees. At knot C 
the jet exhibits an apparent 70-degree kink, and both 
the radio and optical apparent MFPA rotate to remain 
broadly parallel to the new jet direction. The apparent 
MFPA differences in the centre of knot C are consistent 
with zero, although there is a small region to the south- 
west of the knot where the differences increase to ~ 25 de- 
grees and the errors are 8 degrees. On the north-east side 
of knot C, the differences are about —18 ± 11 degrees. In 
knot D, the misalignment increases to —24 ± 11 degrees. 
The asymmetric 'ridge' of radio polarization mentioned 
previously (to the north-east of knot D) does not have 
a simple polarization structure, although here too the 
MFPA is broadly parallel to the jet. 



3.3 Motivation for an oblique shock model 

iKirk fc Heavens! (jl989r ) showed that the spectral index of 
radiation from an oblique shock is flatter than that ex- 
pected from a parallel shock, which is consistent with our 
observation of the radio spectral index an in the range 
0.4-0.6 at knot C. However, the strongest evidence for 
the oblique shock can be found in the polarization maps 
of the jet. The available optical and radio polarimetry 
(Figure f2]) indicates a highly-aligned magnetic field at 
knot C, whe re the jet exhibits a la r ge ap parent change 
of direction. IWorrall fc Birkinshawl (|2005f ) outlined how 
such a change could arise if an oblique shock was present 
at that point, and suggested that such a shock could be 
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Figure 2. Top Panel (a): Optical intensity contours and polarization vectors (rotated through 90 degrees). Contour levels 
are 0.08 X (1, 2, 4, 8, 16, 32, 64) electrons sec -1 , and a vector of length 0.15 arcsec corresponds to 100 per cent polarization. 
The colour scale shows the ratio between percentage of polarized flux detected in the radio and optical bands, where zero 
indicates relatively little polarized radio flux. Bottom Panel (b): Radio intensity contours and (rotated) polarization vectors 
to the same angular scale as (a). Contour levels are 0.36 X (1, 2, 4, 8, 16, 32) mjy beam -1 , and a vector of length 0.15 arcsec 
corresponds to 100 per cent polarization. The colour scale shows the difference between the optical and radio magnetic field 
position angles in degrees. 
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Figure 3. A three-dimensional perspective view of the ge- 
ometry used to d escribe the oblique shock in 3C 346, after 
ICawthornd (2006). The solid grey area is the plane of the sky 
for an observer viewing along the —k direction, and the plane 
of the shock contains the origin. The jet (indicated by the 
shaded cylinder) is directed towards the observer, where the 
upstream and downstream jet components are split by the 
shock front. See text for a description of the labelled angles. 



caused by the turbule nt wake of the companion ga laxy to 
3C 346 (see figure 9 of lWorrall fc Birk inshaw 2005). Con- 
sidering the proposed angle of the shock, they predicted 
a relatively simple polarization pattern at this point. Us- 
ing the new polarization measurements, we now go on to 
explore the properties of such a model and to examine 
how well it fits our observations. 



The geometry is indicated in Figure [3j where the 
shock is shown at the origin. The upstream flow direc- 
tion is shown along the — z-axis, and the plane of the 
shock front is inclined at an angle r\ to the upstream flow. 
The bulk velocity of the upstream plasma is /3 u c, and the 
downstream velocity f3d C and actual deflection a ngle of 
the flow £ were given bv lLind fc Blandfordl l|l985h as: 



Pd = 



and 



[(1 



/3 2 cos 2 77) 2 + 9/3* cos 2 rj sin 2 r\ 



1/2 



tan£ 



tan 2 V (3(ll 



3/3„ sin r) 



(1 



(1) 



(2) 



tan7/(tan 2 rj + 1 + 2/32) 

The plane formed by the upstream and downstream 
flow vectors is rotated by a counterclockwise angle e rela- 
tive to the :rz-plane (when viewed from the +z-direction). 
This structure is viewed by an observer in the plane of 
the sky along the —k direction, which is also in the xz- 
plane and which makes an angle 9 with the 2-axis. The 
apparent deflection angle of the flow in the plane of the 
sky is ip, which is given by 

— sin £ sin e 



tan ip = 



(3) 



— sgn(sin 9) (sin £ cos e cos 9 — cos £ sin 9) 

where the arctangent should be evaluated using all four 
quadrants, and the angle ip as defined above is the coun- 
terclockwise angle from the projected z-axis towards the 
projected downstream flow direction. Figure [3] shows ip, 
the 360-degree complement to ip. 

I n the following s ections we examine two models 
from ICawthornd (|2006l ') and apply them to the 3C 346 
data to see how the magnetic field geometry may give 
rise to the observed polarization at the shock. We use 
the apparent deflection angle ip as a constraint to exam- 
ine the distribution of the free parameters in the model: 
the obliquity of the shock (angle 77), the viewing angles e 
and 9, and the upstream flow speed /3 U . 



4 MODELLING THE SHOCK 

Since the available data on the 3C 346 jet strongly sug- 
gests the presence of an oblique shock at knot C, we have 
examined this region using models of the polarization 
from o blique and conical shocks pres ented bvlCawthorne 
(2006). The se models were based onlCawthorne fc Cobb 
(1990J) and iLind fc Blandfordl (|l985l ). The key assump- 
tions of the models are that: 

(i) the shock front is a stationary surface; 

(ii) the equation of state (T = 4/3) is that of a rela- 
tivistic gas, and the magnetic field is dynamically unim- 
portant; 

(iii) the emission is from a thin layer of material just 
downstream of the shock, so that radiation transfer may 
be neglected; 

(iv) the shock structure does not evolve over at least 
one light-crossing time; 

(v) the emitting material is optically thin; 

(vi) any emission from upstream material is negligible 
compared to the increase in synchrotron emissivity from 
material at the shock front. 



4.1 Single-component field 

The first model (presented bv lCawthorne fc Cobb|[l990l ) 
assumes that there is no ordered magnetic field present 
upstream of the shock (i.e. any upstream field is com- 
pletely tangled and any emissions from upstream mate- 
rial are unpolarized) . The polarization at the shock itself 
would then arise entirely from the tangled upstream field 
becoming partially ordered due to compression by the 
shock. In the 3C 346 jet, significant polarization is ob- 
served upstream of knot C, but we have nevertheless ex- 
amined this model to see how well it might reproduce 
the observed downstream polarization properties. The 
polarization angle \t of synchrotron radiation from the 
shock-compressed tangled-field component was given by 
ICawthornel (|2006h as 

cos e cos ri(Pu — cos 9) — sin n sin 9 , . 

tan Xt = ri { TrT 4 

sinecos 77(1 — p u cos 9) 

where \t ls the electric- vector position angle (EVPA), 
measured counterclockwise from the j/-axis in the plane 
of the sky. The fra ctional polarizatio n of this shocked- 
tangled component (|Cawthornell2006l ) is 
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16 43 48.80 48.75 48.70 48.65 48.60 

RIGHT ASCENSION (J2000) 



Figure 4. The geometry of the single-component model as applied to the 3C 346 
data, with axes and polarization vectors superposed. The electric-vector position 
angle is measured with respect to the yy-axis and is denoted by \. The electric vector 
polarization E is indicated by the dot-dash orange line, and the magnetic vector B 
by the solid orange line. This sketch should be used as a guideline only. 



m t — mo- 



(1- 



(5) 



where mo = (a + l)/(a + 5/3) is the fractional polar- 
ization from a uniform field, which depends on the spec- 
tral index a of optically thin synchrotron emission. The 
compression coefficient (k, the ratio of the upstream to 
downstream number d ensities) for an upst r eam L orentz 
factor 7„ was given bv lCawthorne fc Cobbl (| 199dh as 



1/2 



~? U /3 U sin 77(8/33 sin-* 77 - 7^) 1/2 ' (6) 

The quantity is the angle between the shock nor- 
mal and the line-of-sight in the rest frame of the down- 
stream plasma, and is found by resolving fc into compo- 
nents parallel and perpendicular to the upstream flow 
in its own rest frame (k'» and fc" respectively, from 
ICawthornei 1200^ : using these, the component normal 
to the shock (k 1 ^) can be found and transformed into 
the rest frame of the downstream plasma to give an ex- 
pression for -;/)'. The component normal to the shock is 
k'n = k'a sin rj" — fc" cos rj" , where tan 77" — 7„ tan 77, fcy' = 
(cos 6" — /3 u)/(i — PuCos9) a nd fc" = sin#cose/(7 u [l — 
Pu COS 6»]) (|Cawthornel l2006h . In the rest frame of the 
upstream flow, the downstream velocity is expressed as 
P' d ' = SPu sin?7"/2 - l/(2p u sin 77"), and so finally 

k n + Pd 



cosip' — 



1 + K/3" ■ 



(7) 



Tangled Field Component 








X 














\ Pd = 0.962 






f = 5.7° 


6 = 5.5° 


\ 


ij, = -65.4° 


e = 53.1° 


\ 


X* = -59.0° 


■q = 13.7° 




k = 0.42 


P u = 0.980 


\ 


77ii = 0.24 



Figure 5. An example of the solutions allowed using the 
single-component model with only tangled field. The figure 
shows the model viewed in the plane of the sky, with k point- 
ing out of the plane of the paper (cf. Figure [3j. Although this 
model can reproduce our observations (Figure |4j, the geome- 
try of these solutions at such a small angle to the line-of-sight 
is implausible: see text for discussion. 



4.2 Results using the single-component model 

We search for sets of (p u , sin 9, sin 77, cose) parame- 
ters that can reproduce the jet kink and polarization. 
From the radio data, the measured spectral index a is 
0.55 ± 0.05 at the shock (knot C; see Section [3T| . and 
the apparent deflection angle of the radio jet (as defined 
in Section 2| is —67 ± 1 degrees. The observed EVPA at 
the shock is — 57±2 degrees, and the polarization fraction 
m t = 0.25 ± 0.02 (see Figure [4}. With these constraints, 
we step in p u using a step-size of 0.01 from 0.60-0.99, 
and we map the parameter space with a resolution of 



0.005 in the range ^ sin 9 ^ 1, ^ sin 77 ^ 1 and 
— 1 ^ cose ^ 1. Sets of (p u , sin 9, sin 77, cose) parameters 
that can reproduce the observed polarization features are 
stored as possible solutions - otherwise, they are rejected. 

The only solutions allowed by our data for this model 
lie at small angles to the line-of-sight (9 < 6 degrees), 
where the obliquity of the shock is large (77 < 15 de- 
grees) and the upstream plasma is highly relativistic 
(Pu > 0.98). The allowed parameter ranges are given in 
the first section of Table [3] and an example is shown in 
Figure [5] Although the essential features of the shock can 
be reproduced using this model, these solutions lead to 
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Table 3. Th e results of model ling the jet in 3C 346 using 
the models of ICawthornel d2006h . The minimum, best-fit and 
maximum values are shown for each parameter. The best-fit 
values and errors were estimated using the median, lower quar- 
tile and upper quartile of the parameter distribution. See text 
for discussion. 



Parameter Min Best-fit value 



Max 



Single-component model: Tangled upstream 


field 


Pu 


0.98 


99+ 00 


0.99 


sin 9 


0.055 


065+ 025 
U.UDO_ Q00 


0.100 


sin n 


0.150 




0.245 


cos e 


0.570 


605+ 015 
u.ouo_ 015 


0.675 


Two-component model: Uniform and tangled upstream field 




0.60 


91+ 05 

-0.07 


0.99 


sin 9 


0.005 


oA.r\+ - 135 
u.z^u_ 115 


0.655 


sin n 


0.105 


+0.105 
U -' '°-0.135 


0.995 


cos e 


-0.660 


0.070ig : «S 


0.995 


/ 


0.21 




1.00 


9 


0.00 


12+ 009 


0.98 


It (mjy) 


1.0 


47.8±£? 


54.3 


Qt (mjy) 


i. 1 


re 1+0.5 
°- i -0.3 


9.8 


Ut (mjy) 


-9.5 


^■ y -2.0 


28.4 


Iu (mjy) 


0.002 


6 5+ 4 - 7 
°-°-2.8 


53.3 


Qu (mjy) 


-5.3 


-0.7i° j 


-0.0002 


U u (mjy) 


-38.0 




-0.001 



implausible geometries: such a small angle to the line-of- 
sight would mean that the nearest upstream radio knot 
(knot B) is at least 15 kpc from the shock. Also, because 
3C 346 is an FR II radio source, the core should ap- 
pear more like a quasar than a radio galaxy, with broad 
emission lines, if the jet follows a straight path at this 
angle all the way back to the nucleus. This is inconsis - 
tent with observations (see | Worrall fc Bi rkinshawl l200ll ). 
IWorrall fc Birkinshawl \20o£ also showed that an angle 
to the line-of-sight 9 ~ 20 degrees was compatible with 
the large-scale bending of the jet due to pressure gradi- 
ents in the cluster gas. Since we believe the cluster atmo- 
sphere produces the large-scale gradual bending of the 
flow from the core out to knot E, we reject the single- 
component model and instead investigate improved mod- 
els which take into account the upstream magnetic field. 



4.3 Two-component field 

The improved model treats the magnetic field upstream 
of the shock as two components. As in Section 14.11 one 
component of the upstream field is completely tangled, 
and contributes to the observed polarization only after 
becoming ordered at the shock-front. The remainder of 
the magnetic energy is present in a uniform field parallel 



to the upstream flow: since this field is frozen into the 
plasma, the uniform-field component will remain paral lel 
to the jet in the downstream region <|Cawthornell2006r i. 

The electric vector polarization position angle of the 
uniform- field component, \ u , at the shock is given by 



tanxu 



sin £ sin e 



sin £ cos e cos 9 — cos £ sin 9 ' 



(8) 



and since the field is uniform, the fractional polarization 
of this component will simply be mo. 

We examine the distributions of Stokes parameters 
for the shocked-tangied-component (It, Qt and Ut) and 
uniform- field components (/„, Q u and U u ) in addition 
to the distributions of geometric parameters (9, n, e) and 
the speed of the upstream flow (f3 u ). Using this model, we 
assume that all the emission from the bright knot orig- 
inates from a small region immediately downstream of 
the shock front. In this post-shock region, the observable 
Stokes parameters are denoted I , Q and U (Stokes V 
is zero), and the downstream spectral index is a . The 
Stokes parameters are related by 



(9) 



The intensity due to the shocked uniform-field com- 
ponent is given by 

In = gh, (10) 

where the ratio g is derived in Section 14.3.11 Since we 
infer the value of \u from the jet direction (Equation [8}, 
we then use 

I u rn 



Io \ 




t Iu " 




Qo 


H 


Qu 




Uo J 




v u u , 





}u = ±- 



y/l + tan 2 2 X u 
and 

U u = Qu tan 2\u 



(11) 



(12) 



to find the polarization from the uniform-field compo- 
nent in the post-shock region. The polarization from 
the shocked-tangled component was then obtained using 
Equation [9] 



4-3.1 The fractional intensity of the uniform-field 
component 

The fraction of Stokes / from the uniform-field compo- 
nent in the downstream (post-shock) region is defined 
using 

Iu Iu 
g= T = I u +It' 



(13) 



Using equation 10 of ICawthornel {2006), 
Iu cx (B Pi[ jsinri') 1+ao , 



(14) 



where a is the observed post-shock spectral index, B p ^ 
is the downstream magnetic field parallel to the jet, and 
fi' is the angle between the line of sight and magnetic field 
in the rest frame of the emitting gas. The angle between 
the d ownstream flow a nd line of sight in the shock frame 
is f2 (|Cawthornell2006T ). given by 
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cos Q — sin 6 sin £ cos e + cos 9 cos f , 



(15) 



which may be transformed into the rest frame of the 
downstream flow to find fi' using 

cos Q. ~ fid 



cosf2' 



1; cos Q 



(16) 



The magnetic f i eld B v ,d can be found using equa- 
tion 8 of lCawthornel l|200rJ l: 



Bp.d = B p 



1 + ft 2 7„ tan 2 77 
k 2 (1 + H tan 2 rf) 



= hB v 



(17) 



where h has been introduced fo r convenience. Using equa- 
tion 14 from ICawthornd l|2006h . we assume that 



It oc 



^- (2 



(18) 



(with tp' given by Equation [7| will be a reasonable ap- 
proximation for I t for values of a close to 1. In Equa- 
tions [17] and [TF] Bp and Bt are the parallel and tangled 
components of upstream magnetic field, respectively. In- 
troducing / as the fraction of energy in ordered upstream 
field (see Section \4. 3. 2 [) . 



Bp = f [Bp + B 2 ) , 
and 

B 2 t = (1 - /) (-Bp + B\) , 
by definition. After some algebra, we obtain 
1 



9 = 



(19) 
(20) 

(21) 



1 + 



1-/ 2-sin 2 y>'(l-K 2 ) 
/ ' 3k 2 h 2 sin 2 W 



4-3.2 The fraction of energy in ordered magnetic field 



Following similar methods to those in Section 14.3.11 the 
fraction of energy in ordered magnetic field is found using 
the upstream polarization fraction Pi and spectral index 



the upstream magnetic field position angle (MFPA) is 
parallel to the upstream flow (see Figure [T}, which is the 
key assumption that allows us to consider this model. 
The upstream polarization fraction is Pi — 0.22 ± 0.02, 
and spectral index on = 0.7±0.1. Downstream, the spec- 
tral index is a = 0.55 ± 0.05, and the integrated Stokes 
parameters are I = 54 ± 0.091 mjy, Q = 4.4 ± 0.2 mjy 
and U = —9.5 ±0.1 mjy. The apparent deflection angle 
is now tjj = —50 ± 2 degrees, and with respect to the new 
i/-axis, we infer the EVPA of the uniform-field compo- 
nent to be \ u — —49 ± 3 degrees. We therefore choose 
the negative root for Q u in Equation 1111 

As in Section 14.21 sets of (p u , sin 8, sin 7/, cose) 
parameters that can reproduce the observed polariza- 
tion features are stored as possible solutions. The results 
of mapping the parameter space of the two-component 
model within the ranges given in Section 14.21 are shown 
in the second section of Table [3] 

The distribution of (3 U is now very wide, ranging from 
Pu = 0.60 to f3 u = 0.99 - no upper limit to f3 u was found. 
The median value indicates that the upstream flow is 
highly relativistic, with (3 U = 0.9lto o7. The maximum 
allowed value for the angle to the line of sight (9) is 41 
degrees, and while the distribution is not strongly peaked, 
the median value is significantly lower at 6 = 14^* de- 
grees. The distribution of the obliquity angle 77 has a 
sharper peak with a modal value close to 77 ~ 64 de- 
grees, although the median is at 77 — 51 ± 11 degrees. 
The viewing angle e is not well-constrained, since it has 
a distribution with a large spread around e = 86^21 de- 
grees, and a maximum allowed value near 130 degrees. 
The fraction of energy upstream in ordered magnetic 



field is / = 0.319] 



;, and ranges from 0.2 to 1.0, al- 



though the distribution of g (the fractional intensity from 
the uniform-field component) is peaked strongly around 
low values, such that g — 0.1201q'q52- A valid example 
solution that lies within all constraints and errors has 
fl u = 0.92, 6 = 20 degrees, 77 = 55 degrees and e = 95 
degrees. Other solutions exist close to the median values 
given in Table [3] 



/ 



l+sin 2 0' [-L (ipa - l)] TH 



(22) 



where 6' is the angle between the line of sight and the 
magnetic field in the rest frame of the emitting gas, which 
is obtained using 

, cos e - p u 



1 — Pu cos 9 

The parameter s a depends on cti according to 



(23) 



(24) 



Sa ~ 2 r (3 + a) 

where F is the Gamma function. 

4.4 Results using the two-component model 

To apply the two-component model to the oblique shock 
in 3C 346, we treat the upstream flow as travelling to- 
wards knot C from the centre of the X-ray-emitting region 
(knot Xbc), rather than from radio knot B: this ensures 



5 DISCUSSION 

We have shown that the data support the presence of an 
oblique shock in the jet of 3C 346 (see Section I3.3|l . The 
first model we applied to the shock required an unreal- 
istic geome t ry to fit our data. The improved model of 
ICawthornd (|2006T l produced better results, although our 
simulations were not able to provide unambiguous solu- 
tions to the geometry of the shock in 3C 346. Our results 
suggest that the jet (before the shock) is being viewed at 
an angle of roughly 20 degrees to the line of sight, and 
that the jet speed before the shock is highly relativistic, 
with likely values close to p u = O-Qlj^'gij. 

The models described in the preceding sections make 
important assumptions about the state of the shock, as 
listed in Section f?] The first of these, the assumption 
that the shock front is a static surface, means that we 
are unable to use the X-ray emission lifetime to place 
any additional constraints on the shock geometry. Since 
the X-ray emission is localised in a region which is unre- 
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solved using Chandra, it would be possible, in principle, 
to use the X-ray emission lifetime to constrain the geom- 
etry: this would require the X-ray emitting plasma at the 
shock front to be moving at an apparent velocity such 
that the time taken to traverse an angular distance of 
order the size of the Chandra PSF is greater than the 
X-ray emission lifetime in th e rest frame of the observer. 
IWorrall fc B irkinshawl l|2005h reported synchrotron emis- 
sion from the X-ray knot, where the best-fitting X-ray 
lifetime is kS 2 for a given Doppler-beaming factor 8, and 
k = 4.5 years. Assuming a maximum linear size for the 
X-ray emission a, the required condition is 

a_ & sin 0(1 -& 2 ) 

kc " (1-&COS0) 3 ' 1 ' 

where j3 s is the shock speed. However, because the shock 
is assumed to be static (/3„ = 0) , this argument cannot be 
used. We note that the peak of the X-ray emission shows a 
significant offset to the radio and optical peak positions, 
with a separation 0.29 ± 0.06 arcsec (0.80 ± 0.17 kpc). 
Similar offsets have been observed in other kpc-sc ale 
jets (e.g. iDulwich et all 120071 : lHardcastle et al.ll200ll . in 
3C 15 and 3C 66B, respectively), although their under- 
lying cause is currently not well-understood. Even if the 
shock speed was included in the model and constrained 
using the X-ray emission lifetime, the geometry would 
remain under-determined, since there are still more free 
parameters than observables. Another possibility is that 
knot B and C may both represent a complex of knots - 
unresolved at the redshift of 3C 346 - any one of which 
could be the X-ray emitter. 

The other key assumption of the model is that the or- 
dered component of magnetic field upstream of the shock 
is parallel to the upstream flow, and remains parallel to 
the downstream flow after passing through the shock. For 
the MFPA to be parallel to the upstream flow in 3C 346, 
we must assume that /3 U is directed at radio knot C from 
the upstream X-ray knot, rather than from radio knot B. 
It seems reasonable that radio knot C and the X-ray 
emission form part of the same structure, and since the 
X-ray peak falls between knots B and C in our highest- 
resolution radio image, we feel justified in using the X-ray 
knot to define the upstream flow direction relative to the 
shock at knot C. The upstream radio jet also shows sig- 
nificant curvature from the core through knots A1-A3 
and B, so we cannot determine a conclusive /3 U direction 
using the radio data alone. On the downstream side, ra- 
dio knots C, D and E show a good linear alignment, so we 
can determine the direction of (3d with some confidence. 

The MFPA upstream of knot C is also open to in- 
terpretation, depending on whether the optical or radio 
data are used. As shown in Figure the region between 
knots B and C shows the biggest differences in the optical 
and radio polarization position angles, which can be as 
large as ~ 20 degrees. At the highest angular resolution, 
our 0.08-arcsec, 22.5 GHz radio polarimetry shows the 
apparent field to be parallel to the upstream flow in the 
region where the X-ray emission is brightest, although 
the optical data show a significa nt difference in polar- 
ization angle of 10 — 15 degrees. IDulwich et ail (|2007t ) 
showed that such differences may arise due to a strati- 
fied jet, where different components moving at different 



speeds and emitting in different wavebands could give 
rise to apparent differences in the polarization position 
angle of the radiation from each component, even if the 
magnetic field in each is intrinsically the same. A similar 
phenomenon could be occurring upstream of the strong 
shock in 3C 346, if velocity gradients across the jet are 
large and give rise to a spine-sheath structure in this re- 
gion. This is consistent with our observation that the frac- 
tional polarization in the optical is less than that seen in 
the higher-resolution radio data between knots B and C, 
which may result if the optical emission originates on 
smaller scales and is more prone to the effects of beam- 
averaging. 

The polarization position angles at knot C are well- 
aligned in our radio and optical data, so we have used the 
MFPA at this point to determine the apparent magnetic 
field at the shock front. The most likely geometry of the 
shock, as determined by our modelling in Section 14.41 
is in good agreement w i th th e jet geometry adopted 
bvlWorrall fc Birkinshawl (120051). Our resul ts agre e with 
iGiovannini et al.l (|20oif ) and ICotton et al.l (|l995h . who 
argued that the angle to the line of sight (&) should be 
less than about 30 degrees, based o n jet-sidedness and the 
observ ed properties of the core. IWorrall fc Birkinshawl 
(|2005l ) also analysed the galaxy-gas and cluster-gas prop- 
erties using the X-ray radial profile, and observed that 
the pressures of the two gas components would be equal 
at the position of the X-ray bright knot if = 18+f 2 de- 
grees, which is broadly consistent with our median value 
(A = 14±? degrees). 



6 SUMMARY 

We have presented multi-band observations of the jet in 
3C 346 using high-resolution data from the VLA, HST 
and Chandra observatories. The jet has similar optical 
and radio morphologies, while the X-ray emission from 
the shocked region is significantly offset from the optical 
and radio peak positions. The optical and radio polarime- 
try show close correspondence in the apparent magnetic 
field position angle and the fractional polarization at 
knot C, where the jet undergoes a large apparent kink 
of almost 70 degrees in the optical and radio images, and 
which we interpret as an oblique shock. We have explored 
m odels of the polariz ation from oblique shocks presented 
by ICawthornel |2006) to understand the geometry of the 
3C 346 jet, and we justify the key assumptions required 
to apply the models to our data. Although our modelling 
has not provided a unique description of the jet in 3C 346, 
we find that the upstream flow is likely to be moving with 
velocity f3 u = 0.91+0$, where the plane of the shock 
front is inclined at an angle of?? = 51±ll degrees to 
the upstream flow which is at an angle 8 — 14+® degrees 
to our line of sight. It has been necessary to include the 
polarization observed upstream of the shock and model 
the initial magnetic field using tangled- and uniform-field 
components, since models that treat the upstream field 
as completely tangled fail to give a reasonable description 
of the jet. 
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